X 



The thermodynamic hmit for matter 
interacting with Coulomb forces and with the 
quantized electromagnetic field: 



q: I. The lower bound 



Elliott H. Lieb^ and Michael Loss^ 

1. Departments of Mathematics and Physics, Jadwin Hah, 
Princeton University, P. O. Box 708, Princeton, NJ 08544 



' 2. School of Mathematics, Georgia Tech, Atlanta, GA 30332 

o ■ 

g ; July 31, 2004 

O 

I ■ Abstract 

' The proof of the existence of the thermodynamic limit for electrons and nuclei 

interacting via the Coulomb potential, in the framework of non-relativistic quantum 
mechanics, was accomplished decades ago. This result did not take account of inter- 
actions caused by magnetic fields, however, (the spin-spin interaction, in particular) 
or of the quantized nature of the electromagnetic field. Recent progress has made it 
possible to undertake such a proof in the context of non-relativistic QED. This paper 
contains one part of such a proof by giving a lower bound to the free energy which is 
proportional to the number of particles and which takes account of the fact that the 
field, unlike the particles, is never confined to a finite volume. In the earlier proof the 
lower bound was a 'two line' corollary of the 'stability of matter'. In QED the proof is 
much more complicated. 
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1 Introduction 

Some years ago the problem of proving the existence of the thermodynamic hmit for electrons, 
nuclei and other particles interacting via Coulomb forces was settled in the context of the 
non-relativistic Schrodinger equation jTUj. The key ingredients in this proof, in broad outline, 
were: 

a) The stability of matter of the second kind j21 (i.e., a lower bound on the ground state 
energy proportional to the number of particles), which led to an upper bound on the partition 
function Z, and hence a domain independent lower bound on /, the free energy per particle. 

b) A rigorous version of screening together with a variational argument for a lower bound 
on Z, which led to the fact that / could only decrease (with the density p and inverse 
temperature f3 = l/ksT fixed) as the size of the domain Q containing the particles increases. 
Charge neutrality is needed for this monotonicity of / (but not for the lower bound). Since 
/ is bounded, this monotonicity guarantees that / has a limit as \ fl\, the volume of fl, tends 
to infinity. 

Since then much progress has been made in understanding non-relativistic quantum elec- 
trodynamics (QED) and it seems appropriate now to try to extend the proof of the thermo- 
dynamic limit to the QED case. This is not just an idle exercise, for several new matters of a 
physical nature, as well as a mathematical nature, arise. Among these is the fact this model 
completely takes account of everything that we know about low energy physics, except for 
the hyperfine interaction (for which nuclear physics is necessary, as we explain below), and 
except for the fact that the dynamics of the particles (but not the electromagnetic field) 
is non-relativistic. Indeed, no completely satisfactory relativistic Hamiltonian is presently 
available and, therefore, the fully relativistic generalization will have to await further devel- 
opments. Another problem, which is yet to be resolved, is the renormalization of physical 
parameters in order to deal with the infinities that arise as A, the ultraviolet cutoff on the 
electromagnetic field, tends to infinity 

Otherwise, the theory is potentially complete, as we said, and an example of this com- 
pleteness is that it is not necessary to exclude the spin-spin inter-electron magnetic inter- 
action, as in pUj. The usual non-QED approximation is to mimic the interaction by a 
spin-dependent potential, which cannot possibly be stable, and which is, therefore, omitted 
from discussion unless a hard core interaction is introduced to stabilize it. In contrast, in a 
full theory in which the magnetic field B{x) is a dynamical variable and the particles interact 
with the field via a a- B{x) term (but without any explicit spin-spin interaction) is perfectly 
well behaved and stable and has all the right physics in the classical limit. 
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(We note in passing that stability of matter requires more than just the field energy to 
stabilize the a ■ B{x) terms. It also requires the 'kinetic' energy terms {p + eA{x)/c)'^ to 
control the a ■ B{x) terms, and thereby stabilize the system. In other words, the terms 
p-A{x) +A{x)'^ are essential for understanding the interaction of particles with each other at 
small distances; the dipole-dipole approximation while correct at large distances, is certainly 
inadequate at short distances.) 

Another major difference between the Schrodinger and the QED theories of the thermo- 
dynamic limit is the necessity of treating the thermodynamics of the field correctly. In 1900 
Planck ^3] gave us the energy density of the pure electromagnetic field at temperature T, 
which implies that the field cannot be confined to the container Q without invoking artificial 
constraints. As we shall explain in detail later, this requires us first to take a limit in which 
the size of the universe U tends to infinity (after subtracting the enormous pure Planck free 
energy) and afterward to take the limit \Q\ oo. Obviously, the subtraction has to be done 
carefully and that is an exercise in itself. 

In this paper we consider topic a) above — the upper bound on Z or lower bound on / 
(after taking the double limit, of course). We shall reserve topic b) for later. In the previous 
work [TU] the upper bound required only a few lines, as we shall explain below, but our QED 
setting presents significant difficulties that have to be overcome. While the analog of the 
Dyson-Lenard lower bound on the energy is known for this QED case (see fHl 121); it is 
far from sufficient for obtaining the upper bound on Z. 

2 Basic Definitions 

There are N electrons with mass m and charge — e. These are fermions with spin 1/2. 
There are also K nuclei with several kinds of masses Mm (with M > 1800 in nature), 
positive charges Ze and statistics (Bose or Fermi) but, in order to simplify the notation, we 
shall assume only one species with charge Ze and mass Mm. The generalization to many 
species is trivial, the only significant point being that that all the nuclei have a Z-value not 
greater than some fixed number Z. We also assume that the nuclei are point charges, the 
generalization to smeared out nuclei being a trivial generalization. 

The arena in which the particles reside is a large region 1] C M'^, of unspecified shape for 
the present purposes, and volume It is a subset of an even larger domain W, the 'universe' 
which, for simplicity we take to be a cube of side length L. The boundary conditions of the 
EM field on dU is, presumably, of no importance, so we take periodic boundary conditions 
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for simplicity (although it has to be noted that changing the boundary conditions on the dU 
changes the total energy (when the temperature is not zero) by an amount far greater than 
the energy contained in Q). One could dispense with the universe U by confining the EM 
field to the box Q, but this would be questionable physically and we shall not do so here. 
The two limits (i.e., with or without the confinement of the field to Q) would be expected 
to yield the same average energy density in the thermodynamic limit, but we prefer to take 
nothing for granted. 
The Hilbert space is 

Ti. = "Helectron ® "^nuclei ® ^ ■, (2-1) 

where T is the photon Fock space in U and Tieiectron is the antisymmetric tensor product 
A^iL^(fi; C^) appropriate for spin 1/2 fermions. Likewise, Tinuciei is an antisymmetric tensor 
product of Aj^^L^(f2; C^) (for fermions) or a symmetric tensor product of K L'^{Q) spaces 
(for bosons) or a mixture of them in the case of several species. A vector in 7i is a function 
of electron coordinates and spins xi, ...,xn',o'i, ■■■,crN and K nuclear coordinates (and 
possibly spins if they are fermions) Ri, Rr with values in J-", i.e., it is a vector in JF that 
depends on the particle coordinates and spins. 

Units : The physical units we shall employ here are 2mc^ for the energy and Ac/2 for 
the length (where Ac = h/mc is the electron Compton wavelength). The dimensionless fine 
structure constant is a = e^/hc (= 1/137 in nature). The electron charge is then —^/a and 
the nuclear charge is Z^fa. 

The total Hamiltonian is 

H = T + aVc + Hf , (2.2) 

where the three terms are the kinetic energy of the particles, the Coulomb potential energy 
and the quantized field energy, which will be explained in detail presently. 

The partition function is given by the trace Z = Tr exp[— and the pure- field partition 
function is given by Zq = Tt exp[—(3Hf], with (3 = l/Zc^T. We are interested in the free 
energy per unit volume 

/ = -ksT lim lim {logZ - logZo} , (2.3) 

with the understanding that we set = Peiectron ■ |^| and K = pnucieus ■ |^| for some fixed 
densities Peiectron and Pnucieus- We denote them, collectively, simply as p. Charge neutrality 
is not assumed. 

Our goal here is to derive a lower bound to /. We do not claim to prove that the limits 
in (j2.3|) exist. For the present purpose they are interpreted as lim sup instead of lim. 
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We now define the various energies in detail. First, the kinetic energies (in units of 2mc^). 

N ^ K 

rp^r^ei^ ynuc = J2Tf{A) + -^^Tj{-ZA) . (2.4) 

The electron kinetic energy operator for each electron is the Pauli operator 

T^{A) = {a ■ {p + A{x))Y = {p + y^A{x))^ + y^a-B{x) , (2.5) 

which is appropriate for a spin 1/2 fermion in the presence of a magnetic vector potential A{x) 
and magnetic field B{x) = curM(x). The operator p is given (in our units) hj p = — zV. 
The subscript i in Tf{A) in ()2.4j) indicates that this operator acts on the coordinates of 
electron i and the x in ()2.4j) is then Xj. Note that in this model the (^-factor of the electron 
is 2. If it were greater than this we would be in serious trouble because then the stability of 
matter would not hold 5J. (Strictly speaking the result in about > 2 holds only for 
classical fields without UV cutoff. With a cutoff one expects stability of the first kind, i.e., 
a finite ground state energy, but not stability of the second kind, i.e., a lower bound that 
is proportional to the number of particles. Although well known QED calculations say that 
the renormalized, effective (/-factor exceeds 2, QED theory always starts with 2, otherwise 
the theory would not be renormalizable [1^].) 

Since the nuclear charge is ^Z^fa we have —ZA in ()2.4|1 . The kinetic energy operator 
in ()2.4j) omits the a ■ B{x) term, i.e., it is 

T,iA) = ip, + V^A{R,))\ (2.6) 

in which Rj is the coordinate of the j*^ nucleus. A nucleus can have a magnetic moment, 
even if its charge is zero (the neutron) but it often has a (/-factor much larger than 2 (e.g, 
(/ ^ 5.5 for a proton). 

A conventional 'physical argument' might be that since the magnetic moment is inversely 
proportional to the mass, the contribution of the magnetic energy to the total energy is small. 
As mentioned before, however, the inclusion of a dipole-dipole interaction has disastrous 
consequences, no matter how small the coupling is. It follows that to include this magnetic 
interaction, and hence to include the hyperfine interaction, we would have to take explicit 
account of the nuclear magnetic form factor. In essence this really means thinking of the 
nuclei as a compound system of quarks. The interaction of quarks with the EM field will 
surely have diamagnetic components. The effects of the nuclear structure, which would add 
a self-energy proportional to K, could be included, but we prefer to avoid this discussion 
here. 



ELML/July 31, 2004 



6 



We use the Coulomb gauge to describe the EM field and its energy. In this gauge only 
the magnetic field is a dynamical variable, i.e., the curl- free part of the electric field is not 
an independent dynamical variable, for it is determined by the particle coordinates and 
Coulomb's law. This choice of gauge is essential because, as we have said elsewhere |S1 ITT|. 
it is the only gauge in which the correct physical EM interactions (including the spin-spin 
magnetic interaction) can be obtained by a minimization principle. 

The field energy (in our units) is 



3 

^/ = EEl^l«K^K(^) ' (2-7) 

k j=l 

where the three operators Si{k) = {ai{k), 02 (/c), 03 (/c)) are boson annihilation operators of 
momentum k satisfying the canonical commutation relations [aj(g) , a* (A;)] = Sij6g^k, etc. 
(We are using the convention of three-component quantized fields introduced in ^1], which 
means that we have to subtract the Planck background energy in ()2.3|) for three modes 
instead of two. The advantage of this formalism is that we do not have to introduce the 
fc-dependent classical polarization vectors ex{k) to insure the 'divergence-free condition' on 
the vector potential A{x).) The k sum in fj2.7|) is over /c's of the form k = {2tt/ L){ni, n2, n^) 
with integer n^, but A; = is excluded. As L — * 00, (27r/L)^ times the sum over k becomes 
an integral Jj^g, but for a finite L we have to do the sum carefully in order to get the correct 
cancellation. 

The vector potential is obtained by first defining the vector field 

C{x) = y[t) 5Z ^a(^)^ - a*(fc)e-*^] , (2.8) 

^ \ J ^ II 

where Xa ^ 1 is a radial function that vanishes outside a ball of radius A. Then 

A{x) = i curl C(x) . 

In the Coulomb gauge the electrostatic energy aVc is given by a simple coordinate- 
dependent potential 

It— t I~'~^ I/? — f?l ^ It- — Rl ' (2-9) 

l<i<j<N ' * ^' l<i<j<K ' ^' l<i<Nl<j<K' * 

Having introduced the free energy and its component parts we can now discuss the 
physical and mathematical problems addressed in this paper, namely the difference between 
the QED problem and the non-QED problem in [TUj . 
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The non-QED problem: In this case there is no need to introduce the 'universe' 
U, the field energy Hf or to take the double limit in ()2.3|1 because there is no particle-field 
interaction via the field A{x). To obtain a lower bound to / one simply writes ^10, Theorem 
2.2] H = H1 + H2 with 

Hi = ^T + aV, ^2 = ^T. (2.10) 

One then bounds Hi from below by —c{N + K), where c is a universal constant, which 
follows from the stability of matter bound. (The constant c changes when we replace T by 
T/2 but it is always finite.) Then, we can bound / by 

/ > -c(Pelectron + Pnucleus) " fc^T lim log Tr e"'^'^/^ . (2.11) 

s2— >oo 

The latter trace is just the partition function of an ideal gas of twice the mass and has, as 
the elementary textbooks tell us, a finite / (which can be easily bounded from below by a 
shape independent /). Thus, a satisfactory lower bound to / is a 'two-line argument' in this 
case. 

The QED problem: Several points have to be considered. 

1. The precise cancellation of the background Planck field energy, which is enormously 
greater than /, has to be done carefully. 

2. We know from earlier work TB] that the stability of matter does not hold for the Pauli 
kinetic energy unless the field energy is added to the Hamiltonian. This implies that 
we must somehow borrow some field energy to stabilize the Pauli analog of Hi above, 
but not too much to spoil the delicate cancellation in item^ The situation with the 
quantized field energy Hf in place of a classical field energy f^^ B{x)'^dx/8Tr is even 
more delicate; this extension was first made in |2| using the results in ^Hl- The idea 
in [2] is to replace the field energy term by one that is localized near the nuclei. Our 
approach will be somewhat different and will only involve localization of the electrons, 
in a manner reminiscent of the original Dyson-Lenard proof [3 of stability of matter. 
(However, if we do not care that our lower bound for / does not go to zero as p — > 
then localization is not needed - see sectional) In any event, the stability of matter 
using the Pauli operator and Hf requires a bound on a and Zo? , when there is no 
ultraviolet cutoff; this is a feature not encountered in the non-QED problem. We shall 
prove that / is bounded below (for all A if a < 1/98 and if ZcP' < 1/468 (corresponding 
to Z < 39 for a = 1/137. These bounds can certainly be improved, with some effort, 
but we shall not attempt to do so. 
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3. The lower bound on / should be shown to have a classical limit as the nuclear mass M 
tends to infinity, independent of the statistics of the nuclei. 'Classical' means that the 
dependence of / on M has the form / ~ Kp^^ciensksT [| log(/5/M) + logpnucicus - l] • 
This would coincide with the experimental observation that the nuclei may as well be 
considered to be fixed in space. 

4. While an infrared cutoff is not needed for our considerations, the ultraviolet cutoff A 
is essential. One would like to take the limit A ^ oo after a renormalization of the 
electron mass and possibly its charge. At present, it is not known how to carry out 
this program, although some primitive steps were taken in [T2j. In any case, in order 
to show that A ^ oo limit can be taken for the thermodynamics it is appropriate and 
necessary to show the following. There are three functions fadP, p, Z), (/nucleus (A, 2') 
and 5'eiectron(A) such that the free energy per unit volume, p, Z, A), which depends 
on all four parameters, can be decomposed as 

f{f3,p,Z,A) = U{f3,p,Z) + 5'clcctron(A) Pelectron ~l" 5'nucleus(A, Z^ Pnucleus • 

(2.12) 

We make the additional requirement that /oo(/5, p, Z) oc p as p 0. With no particles 
there should be no free energy. (It is also physically desirable that the self-energy term 
(?eiectron(A) should bc large and positive when A is large. We succeed only partially in 
this respect, as discussed in the remark following Theorem 13.11 ) 

In this paper a lower bound of the form (|2.12p will be derived and it will have the indicated 
properties. This is done is in section 01 with help from section 0] and the Appendices 

If we forego the property that /oo(/5,p, Z) oc p as p ^ then it is significantly easier to 
obtain a lower bound of the form ()2.12|) . This is shown in sectional 

3 Main Theorem and Proof 

In this section we show how to decompose if in a manner reminiscent of (j2.1Up and how 
to use this to prove a lower bound of the form ()2.12|) . under the assumption that certain 
inequalities hold. These inequalities will be proved in subsequent sections and appendices. 
In the following ^nucleus is the partition function of non-interacting bosons of mass Mm. As 
M oo this partition function tends to the classical partition function if T > 0. This was 
one of our requirements in item 01 in section |21 

Thus, we shall prove the following in this section. 
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40 / V 468 > 



THEOREM 3.1 (Lower bound on f). Assume that a < 1/98 and Za^ < (§)( 
(This means Z < 39 when a = 1/137. j Then ^2.11]) holds with 



/oo(/3, P, Z) > —kT {peiectron 111 2 + Pelcctron lll(Pelectron) ~ ^electron} 

— /^Tpelectron Ih (S7l(kT)^C^ ^) + /nucleus Pnucleus) ; (3-l> 



<7eiectron(A) > " (128 + 0.0032) A^/^ - (^-1-^ - (|^) (468) (40) A^/^, (3.2) 

5'nucleus(A, ^) > , (3.3) 

where /nucleus is the free energy per unit volume for non-intercting bosonic nuclei of mass 
mM and 

^sa = ^ - «^ max{64.5, ^Z} - ^ exp | ^3^^ } (3.4) 



See equation ()3.39|) for a more general expression. 

REMARK: The bounds ()3.3|) on 5'eiectron(A) and (y'nucieus(A, Z) should actually be large 
and positive owing to the self-energy of the charged particles. In we showed that the 
dependence of this self-energy on A is somewhere between Is?/'^ and A^^/^ and we conjectured 
that A^^/^ was the correct dependence. If A^^/'' is, indeed, correct then we could easily invoke 
the methods of and add a term proportional to A^^/'' to (7eiectron(A), which would dominate 
the A^/^ term and thereby leave us with the desired large positive value for 5'eiectron(A). On 
the other hand, if A'^/^ is correct then adding A'^/^ to 5'eiectron(A) will not produce a positive 
self-energy (for large A, of course). Since the A^/^ - A^^/^ question is not yet settled, we shall 
not burden this paper further with the question of the sign of the self-energy. We note in 
passing, however, that the simpler theorem in section has a huge negative A dependence 
(i.e., — A^), which cannot be compensated even by A^^/''. 

Proof. The first step is to localize the electrons. Decompose Vt into disjoint cubes of side- 
length A^*. There will be approximately \VL\ A^'' of these cubes. For each cube 7^ we take 
a smooth function (pi centered on 7^, whose support is in a cube Vi of twice the side length 
and in such a way that 

^ <pe,{xf = 1 for all a; e 1^ . (3.5) 
We can and do require that |V0£(a;)| < lA"*. 
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It is easily seen that for each point x there can be at most 8 distinct values of i for which 
(f)e{x) 7^ 0, and hence, by the standard IMS localization formula, 

T^(A) = J2M^)T''{A)M^) - J2 |V0£(a;)P > </),(x)T^(A)</),(a;) - 128A2^ . (3.6) 

e i e 

for all A{x). By applying this to each term in yiectron^ recalling ()3.5p . we obtain the 
inequahty 

j^electron > ^ ^^yelectron^^ _ 12SNA^-^ ^ (3.7) 
/ 

where / = En) is a multi-index and = xn) = 4>eiixi) ■ ■ ■ (Pij^^xn) is a 

product function, which satisfies Yli ^AK)"^ = 1 when X e . Here, X = (xi, X2, ■ ■ ■ , xn) 
collectively denotes the N electron coordinates; similarly, R = (Ri, R2, ■ ■ ■ , Rr) denotes the 
K coordinates of the nuclei. 

Armed with this localization, and recalling ()3.5p . we can write (for an arbitrary constant 
Ci>0) 

N 

^electron > {T[' {A) - + V^A(x.) | } $,(X) 

/ i=l 

AT 

+ $^<^'/(^)$^C'ib, + A(x,)|<l>,(X)-128iVA2^ , (3.8) 
and hence our Hamiltonian is bounded below as 

H>H2 :=J2 '^i(^) [H''''^ + - l^SNA^' , (3.9) 

where 

N 

^rad ^ ^ I^P(^) _ ^^1^^ ^ y^A(Xi)|} + (3.10) 
i=l 

TV 

^rel _ ^nucleus + ^ |pi + V^A{x,)\ + aV, . (3.11) 

1=1 

The superscript 'rel' is meant to suggest that (|3.11|) is a Hamiltonian of relativistic-like 
electrons and non-relativistic nuclei. Note that $/(X)ifj$/(X) = Hf since Hf does not 
depend on the electron coordinates X. 

The reason for adding and subtracting the 'relativistic' kinetic energy operator \p + 
y/aA{x)\ is that Tr exp{— T^(74)} is not bounded independent of the vector field A, but 
Trexp{ — |p+ A/ay4(a;)|} is uniformly bounded. 
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In sect. EJ Lemma [4. H we shall show that the sum appearing in ()H.1()|) is bounded below 

by 

-^V^ E ' / B{yrWj{y)dy - NCf , (3.12) 
3v3 J Jm.3 

where Wj, given in ()4.2j) . is a sum of characteristic functions of subsets of Q whose total 
volume is at most 8A^A~'^*. This function has the following properties for each J: 

/ Wi{y)dy<8NA-^% Wj{y)<8. (3.13) 

The constant C stems from the application of a Lieb-Thirring inequality and it is bounded 
above by 0.06003. 

If D{y) is the (vector) operator obtained from the a part of (j2.8|) . namely 

3/2 



(9 \ ' 1 
T) E^a(^)^^(^)^^''' (3.14) 

= ^[t) E ^a(^)^^^''' [^'^(^) - ■ ' (3-15) 



fc 

then Bill) = D{y) + D*{y) and Schwarz's inequality leads to 

B{yf < AD\y)D{y) + 2 [Z}(|/), Z}(|/)*] , (3.16) 

which then implies that 

B{yfWi{y)dy <4 [ D*{y)- D{y)Wj{y)dy + -A^-^'N . (3.17) 

Altogether, (j3.12j) and the definition of H'^^'^ and H^'^^ lead to the lower bound 



H > H2>H3:=^<^i{X) [Hl^'^ + H''^] <l>j(X) - 



128A2^ + Cl + 1^!^a4-3. 

3v37r 



(3.18) 



where Hl'^'^ is a replacement for if ""^^i given by 

= Hf- J2 '^^(^)' / • (3.19) 

3v3 ^ JR3 

/n summary, our lower bound Hamiltonian contains three parts: A constant propor- 
tional to N , a perturbed field energy ^jH^'^^^i, and the Hamiltonian of 
'relativistic' electrons and non-relativistic nuclei. Their definitions depend on a constant 
Ci > and on the parameter s, which defines the electron localization. These will be chosen 
later. 

Our final goal is to prove the following upper bound on Tr exp{— /^i^s}, which will then 
complete the proof of Theorem 13.11 
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LEMMA 3.2. 

1 N 



e-/3C3bl^p 



X 



exp|/3iV(128A2^ + C2 + i^^^£^A4-3« + l^|ln(l-£)|A3-3^)l , (3.20) 
t 3v37r 3 J 

where ^nucleus (-^) is the partition function of K non-interacting bosons of mass M (even if 
some of the nuclei happen to be fermions). The numbers (assumed to be < 1, see iC.10\} } 

IS 

4096 71 

e = ^^C,C. (3.21) 

Proof. First, we dispose of the localization function $(X) that appear in the Hamiltonian 
of (j3.18p . By Lemma IB. II in Appendix B, and the fact that ^i{X_y = I = identity 
operator on our Hilbert space ()2.1|) . 

Tr e-^^^ < Tr e-/^(^r'^+^'") exp l(3N{128A'^ + Cl + 1^!^£i^a4-3.) 1 _ (3 32) 

t 3v37r J 

Next, we introduce another constant < C3 < Ci, and write 

N N 
^rel ^ ^nucleus + ^3 \P^ + V^A{x,) | + J](Ci - C3) \pi + V^A{Xi) \ + aV^ . (3.23) 



i=l 1=1 



Using the result of pr7|, as stated in ^Hj (see also UM,), the last two terms, taken together, are 
positive as an operator on the tensor product space of the nuclei and the spin 1/2 electrons 
provided that 



Ci-C3> max|^— , nZ}a . (3.24) 
^- ^0.032' ^ ^ ' 



If is true then 



N 



^rel > ^rel _ ynucleu. + C'3 ^ + V^A{Xi) \ . (3.25) 



i=l 



Using this and the Golden- Thompson inequality 

Tre-^(^f"+^") < Tre-^(^^"+^") < Tr (^e-^^^'e"^^") . (3.26) 

We evaluate the trace in the Schrodinger representation in which the field A is a c- 
number field and the trace is just integration over this classical field. This is a rigorous 
technique in quantum field theory and we explain it in some detail in Appendix El We 
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then use the fact that for a fixed classical field A{y) the operator e ^^i''^ has a kernel 
exp{—PHl^^}(^,R ; 2^,R), and we can write the right side of (j3.26|) as 

Tre-^^J'e-'^^^" = 2^ j V{A) j dXdR 

e-^^l''(^){X,R,_; X,R,)_ (A|e-^^^^' \A){X) , (3.27) 

where the factor 2^ comes from the electron spins. (Note that has no electron-spin 
dependence.) There might be a factor for nuclear spins, but we shall ignore this subtlety in 
order to keep the notation simple. In any case it can be absorbed in the factor 2nucicus(-^)- 
Obviously, the matrix element e~'^^i'''''"^^(X, R ; X_^K} is the product of a factor depending 
on the nuclear coordinates R and a factor depending on the electron coordinates X. Each of 
these factors is an (A-dependent) partition function of non-interacting particles. The former 
depends on the statistics of the nuclei. If the nuclei are bosons e~^'^i°'*^'^^(X, ^ ; X_,K) can 
be written as 

-^per (^x^{--j-^{p - Z^/a A{x)f}{Rj, Rk)^ ^ M {eyi^{-l3Cz\v + \/aA{x)\ }{xi,Xj)) , 

(3.28) 

where per indicates the permanent of the K x K matrix (indexed by j. A; G {1, . . . , K}) and 
det the determinant of the N x N matrix (indexed by i,j G {1, . . . , N}). If the nuclei are 
also fermions or a mixture of fermions and bosons we have to replace the permanent above 
by the corresponding symmetrized or antisjmametrized product, i.e., permanent for a boson 
species or determinant for a fermion species. 

For our purpose here, namely an upper bound, we may assume from now on that all the 
nuclei are bosons. The reason is that the K x K matrix above is positive definite, and so 
is the N X N matrix. It is a fact that the determinant of a positive definite matrix is not 
greater than the permanent. Indeed, the determinant is less than or equal to the product of 
the diagonal entries while the permanent is greater than or equal to the same product. (See 
(0)-) Since (A|e~^^i'"*|A) is positive, we can use this upper bound on the determinant to 
obtain the following upper bound to the right side of ()3.27p 



j V{A) j dXdR-^pei (e-M(P-^v^^(-))'(/?^.,i?, 

1.-Q (^,-/3C3b+v^^(.)|(^^^^^)j X (A|e-^^^" \A)iX, ) , (3.29) 



i=l 



Since all the factors in (j3.29j) are positive, we can appeal to the diamagnetic inequality 
and delete the field A from the second factor. (The diamagnetic inequality is well known 
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and states that exp[—(3\p+^/aA{x)\'^]{xi,Xj) < exp [— (x^, Xj), which follows from 
the Feynman-Kac representation; it is also true for exp[—f3\p + ^/aA{x)\]{xi,Xj), thanks 
to the fact that e~'^' = e~^~^^^'^^dt/y/7ii.) Note that it was first necessary to replace 
the determinant by the product of its diagonal elements and then to use the diamagnetic 
inequality; otherwise we would have to worry about the minus signs in the determinant 

Similarly, we can set A = in the first (permanent) factor. The reason is that the 
permanent can only increase if we replace each matrix element by its absolute value and 
then replace that, in turn, by a larger number. But the diamagnetic inequality (actually, the 
Wiener integral representation, to be precise) tells us that this is achieved by setting A = 
(even if Rj ^ Rk). 

In this manner we obtain the upper bound 

/"I 1 ^ 

/ dXdR— perm (e-^^^\R„Rk)) x H ^^)) Tr^e~''''^^"(^' ) • (3-30) 

■ ' i=i 

The expression containing the Fock space trace still depends on the variables X. Using 

flC.14|) from Appendix O we see that 

In Tr^e-^^^^' (X, R) < In Tr^e"^^^ + — |ln(l - e) \ A^-^' , (3.31) 

3 

where e is given in (j3.21|) and where (see (j3.13|) ) 

Giy, X) = J2 ^liKfWiiy) < 8 . (3.32) 

Inequality ()3.31|) is true, as shown in ()C.9|) . ()C.14|) provided the criterion M < eK given 
there is satisfied. Since G{y,2Q ^ 8 ^or all y,X_ this criterion is satisfied with e as in ()3.21|1 . 
and this can be achieved by choosing Ci small enough. 

Recalling that the operators p"^ (associated with the nuclei) and \p\ = (associated 
with the electrons) are Dirichlet Laplacians on the domain f2 we get the upper bound 

1 r ^ 
Tr e-^^^ < Zn^eieus {K) — [Tr e~C^'\ Tr^ e"^^/ 

exp|/?iV(128A2^ + Ci2 + ^^CiA^-=^^ + — |ln(l-£)|A3-3«)l . (3.33) 
I 3v37r 3 J 



This proves Lemma [3.21 □ 



The factor 



Tr e-^^'^ 



N 



can be estimated from above, by the Golden-Thompson 



inequality, as 

1 r /• 1^ 



, ^ exp ( 

Pelectron lll(Pelectron) + Pelectron "I" Pelectron ln(87r/3 C3 ) 

(3.34) 
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To prove Theorem 13.11 we have to consider numerical values for our constants. Let us 
collect together the conditions on them, which are 1)3.211) . ()3.24j] . That is 

(4096) (0.06) vr 



3^3 



-aCi = U9Cia <1, (3.35) 



2 

C3 = Ci — max{ , 7rZ}« = Ci — max{64.5, vrZja > . (3.36) 

.032 

This value of C3 is to be inserted into l3.20t using ()3.34|) — assuming that the two conditions 
on Ci, implied by ()3.35|) and ()3.36p . are satisfied. These two conditions set bounds on a and 
on These are (149) (64.5)^2 = 9613^2 < i and (149)7rZ«2 = (468) Za"^ < 1, as stated 

in Theorem 13.11 

The free constants to be determined are Ci and s. The other constants e and C3 are 
in fl3.21|) and ()3.36|) . respectively. The factor 128aC7i£/3v^7r in ()3.20|) can be replaced by 
0.47/149 = 0.0032 since Cia < 1/149. Our bound is then 

—Pf <Pelectron lu 2 — /? /nucleus (3.37) 

+ /3peiectron (i28A2^ + C^, + 0.0032A^-3^ + ^| ln(l - £)|A3-3^^ (3.38) 

Pelcctron In(pclectron) + Pelcctron + Pelcctron ln(87r/5 Cg ) , (3.39) 

where /nucleus = —kT\Q\~^ In ^nucleus is the free energy per unit volume for non-interacting 
bosonic nuclei of mass mM. 
If we choose, for example, 

, = 149C,« = l-exp(^-Z^51_), (3.40) 

and s = 4/5 in order that the two largest A exponents in ()3.38j) have a common value (8/5), 
and restrict Za'^ < {%){^), then (|niII)-(|S!l) is obtained. □ 



4 Decomposition into Boxes 

In this section we shall give the details of the lower bound, ()3.12|) . of the kinetic energy 
operator contained in ()3.8j) in terms of a Fock space energy operator. We recall the IMS 
localization into disjoint cubes 7^ with side length A~* and overlapping cubes Te with twice 
the side length introduced in ()3.5|) - ()3.7|) . 
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LEMMA 4.1. On the Hilhert space Af^^L'^{^^;C'^) of N electrons with 2 spin states we 
have that for all values of Ci and all vector potentials A{x), 

N 

/ j=i 

-^V^E^^(^)'/ BixrWjix)dx-C!N , (4.1) 

where the function Wi{x), x & Q, is given by 

Wi{x)=Y,Xk{x)<8 , (4.2) 

where Xk is the characteristic function of the cube Tk and where A; G / = (£i, . . . , ^at) means 
that at least one of the ii equals k. Note our convention that each Xk is allowed to enter the 
sum in ( |^.ij| ) at most once, i.e., if k appears 5 times in I then Xk appears once in \4.'/3^ . The 
constant C in \4-l\l is the 7 = 1/2, 3-dimensional Lieb-Thirring constant; C < 0.06003. 

Note: See [7j for the value of C quoted above and see [IHl appendix A] for the fact that 
it is not necessary to include an extra factor of 2 in order to account for the 2 spin states. 

Proof. Fix / = (£i, . . . , ^Tv) and consider the single term 

N 

Yi = [Tf{A) - C^\p, + ^A{x,)\] <^j{X) . (4.3) 

i=i 

In the index set I the index ki appears ni times, the index k2 appears n2 times etc. where 
the numbers nj > 1 and nj = N. 

Our goal is to find a lower bound to (\E', Yj"^) for any in L'^{Q; C^). Let us consider 
the first ni terms in ()4.3p . i.e., {^i"^ ,Yl^Li'^j "^z^); where T is the operator appearing 
in [ ] in ()4.3|) . In evaluating this inner product we can fix the coordinates Xj,aj with 
j = rii + 1, . . . , N and then integrate over them at the end. In other words, the proof of 
our inequality ()4.1|) will follow from the following statement: For each n > 1 and each k, 
every normalized, antisymmetric function ip of n space-spin variables, with support in (F^)" 
satisfies the inequality 



(^ 



n . (4.4) 
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By the arithmetic geometric mean inequahty (a^ + Cf > 2C\\fa) 



>-C^n-C,Ti\2JTf{A)-\p + ^A\ . (4.6) 



Here, denotes the negative part of x (which is always > 0). 
Using the inequahty of Birman and Solomyak P (see also [TS] ) 

- b + v^A|] _ < Tr [ATfiA) - i\p + ^A\Y] '[^ 
= Tr [(3(p + v^A)2 + 4v^a ■ B)]^^^ . 



Tr 



(4.7) 
(4.8) 



By the Lieb-Thirring inequality (but with the added remarks in ^H] to avoid the factor of 
two) ()4.8|1 is bounded above by 

16aC 



3^3 



B(x) dx 



(4.9) 



The bound Wi{X_) < 8 in ()4.2j) comes from the fact that a point x G can lie in at 
most 8 cubes Tf. □ 



5 A Simpler Theorem with a Simpler Proof 

In this section we show how to obtain a lower bound on the free energy per unit volume / 
that is correct in all respects except that it does not vanish as p — ^ 0. Not only is the proof 
simpler but some of the constants are also better. No localization is required. 

THEOREM 5.1 (Simplified lower bound on f). Assume that a < 1/35 and Za^ < 
(||^)(gj-g). (This means Z < 250 when a = 1/137. ^1 Then 

/(/3, p, Z) > —kT {pelectron lu 2 + Pelectron lll(Peloctron) ~ Pclectron} 

- kT Pelectron In {8lT{kTfC;^^) + /nucleus (/?, Pnuclcus) , 



47r(468)(40),3 / 1 



- ( 1 Pelectron " 0.026A^ , (5.1) 

where /nucleus 'is the free energy per unit volume for non-interacting bosonic nuclei of mass 
mM and 

= 18:6 - - 18:6 l-58.5Za^ ^^-'^ 
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/ B{y fdy - NC! . (ISI2I) 

Jn 



See equation for a more general expression. 

Proof. The proof is as in the proof of Theorem 13 . II except that the electrons are not localized 
(but they are confined to the domain Q and the wave function satisfies Dirichlet boundary 
conditions on dQ). In other words, we eliminate $/ and from the equations. The local- 
ization penalty 128 A^A^'^ is eliminated. The function Wi{y) is replaced by the characteristic 
function of the domain Q and ()3.13p is replaced by J^^ W{y)dy = \Q\ and W{y) < 1. 
Expression ()3.12p is replaced, therefore, by 

Thus, we save a factor of 8 because there is no longer a concern about overlapping cells F^. 

The function G{y,2Q is replaced simply by the characteristic function of Q (for all X_), 
whence J^^ G{y,2Qdy = \ fl\. The bound ()3.32|) is replaced by G{y,JC) < 1. In view of this, 
the number e in ()3.2H) and ()C.10|) is reduced by a factor of 8 to 

e' = 87rC4 = 5127iaCiC/3Vs ^M) 

Inequality ()3.17p becomes 

/ BiyYdy <A [ D*iy)- D{y)dy + -A' \n\ . dSTTI) 
Jn Jn ^ 

Lemma l4. II remains true, but with the obvious replacement of ^ J^^ B^Wi by B^ 
in ()4.H) and with ()4.2|) eliminated altogether. 
Lemma 13.21 becomes 



LEMMA 5.2. 

Tre-^^<Z„,eicus(i^)^ 



Tr^e-^^/x 

exp|/?(+ArC? + l^A^|fi| + ||ln(l-.')|A^N)} , (5-3) 

where ^nucicus(-^) is the partition function of K non-interacting bosons of mass M (even if 
some of the nuclei happen to be fermions). The number e' (assumed to be < I) is in 113. 2 
and, as in / l^/. c/6|) . C3 = Ci — max{64.5, vrZja > 0. 

The final task is to choose Ci. Our conditions on e' and on C3 lead, as before, to 
conditions on a and Za"^, namely (18.6)(64.5)a;^ < 1 (or a < 1/35) and (18.6)7rZQ;^ < 1 (or 
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Za^ < 1/58.5). We choose a slightly lower bound for Za'^, namely Za^ < (250/320)(l/58.5), 
and we choose 



We can bound if^^ by 0.47/18.6 = (0.0032)8 = 0.026. These choices lead to Theorem 



A The Schrodinger representation 

In the proof of Lemma (3.21 especially eq. ()3.27p . we evaluated a trace over the full Hilbert 
space in the "Schrodinger representation" in which the field A is regarded as a c-number 
field. For a fuller discussion and justification of this method we can refer, for example, to jHl 
part I, sec. 2], but here we discuss only what is needed in our application. 

First, we note that since the volume of the universe \U\ is finite and there is an ultra- 
violet cutoff A, there are only finitely many photon modes that interact with the electrons 
and nuclei. Each mode is a harmonic oscillator mode and can be described in the usual 
Schrodinger representation by the canonical operators pk and qk, one pair for each /c- value 
and each polarization. In our case the qk is just the Fourier component of A{y) namely A{k). 

The noninteracting modes are infinite in number but they can be ignored since their 
contribution to the trace is easy to compute (Planck's formula). 

In evaluating the trace on the right side of (j3.26|) we can use the q representation for 
the photon modes and the x representation for the space, as usual. The operator Hl"^^ 
involves the electron/nuclei p and x operators but it involves only the qkS and not the pkS. 
On the other hand, H{'^'^ involves the pk and qk operators and the x operators, but it does 
not involve the electron/nuclei p operators. 

Thus, e~'^^i°' is a multiplication operator as far as the photon modes are concerned 
and e~^^^'^'^'^ is a multiplication operator for the electron/nuclei space. In fl3.27p the 
notation dX dR is standard Lebesgue measure while T>{A) means Lebesgue integration over 
the (finitely many) g^'s. It is well known that the (finite) trace of the exponentials of the 
operators that we are considering can be evaluated in the Schrodinger representation by 
taking the x, x and g, q matrix elements and integrating over these variables in this manner. 
Indeed, the kernel {A\e~^^^^'"^\A') is an (X, R dependent) Mehler kernel. 




(5.4) 



□ 
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B A lemma about localization 

LEMMA B.l (localization and convex functions). Let H be a self-adjoint operator 
on some Hilbert spaces Ti with dense domain D{H). Let = 1, ... M be a collection of 
bounded operators such that $j maps D{H) to itself and such that 

M 

5^$*$,=/ (B.l) 

i=i 

Then, for any convex function f{x) and any normalized \I/ G D{H), we have that 

M / M \ 

($,vl>, /(i7)$,vl/) > / 5^ ($,vl>, H^,^) . (B.2) 
j=i \i=i / 

// we also assume that X^^i = ^ ^'^^ ^^^^ f{H) is trace class, then 

M 

TTf{H)>TTf{J2'^*H<^j) . (B.3) 

i=i 

Proof. Note that 



M 



|2 



by applying Jensen's inequality in the spectral representation of H. With Cj = ||$j^|| , we 
have that Xljli ~ II^P — hence we may apply Jensen's inequality once more to 

obtain 



This proves ()B.2|) . To prove ()B.3|) we just sum (|B.2|) over the orthonormal basis {ipn} of 
eigenfunctions of X^jli and obtain 

M M M 

Tr/( J] $*//$,) = J]/((^„,5^$*i/<|.„^„)) = J]/(5^(<|.,^„,i7$,^„)) 

j=l n j=l n j=l 

M M M 

^ E E^*^^-^- /(^)'^^-^-) = Tr 5^ <f */(^)% = Tr 5^ /(i/)$,$* = Tr/(if) . (B.6) 

n j=l j=l j=l 

□ 
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C Perturbed black-body radiation 

One of the problems in section El is to estimate the partition function of the Hamiltonian 
Hl^'^ (see fl3.19|) ) in terms of the partition function of the universe, h(, with an error term 
that depends only on the number of electrons and nuclei and, possibly, on the ultraviolet 
cutoff. 

In general, let us consider a Hamiltonian on Fock space of the form 

Hm = hcY, {h,k'5'''\k\ - M^'l,} a*{k)a,{k) , (C.l) 

k,k' 

with M self adjoint, i.e. M^'-^, = Ml',\. Recall that k = ^(ni, n2, ris) with rii an integer (but 
k = is omitted) and i,j G {1, 2, 3, }. The volume of U is L^. 

The matrix {•} appearing in (jdlj) . can be written, in an obvious notation, as K — M. 

Later on we shall present the M we are interested in for the purposes of this paper, but 
for the moment let us consider the partition function 

In^M = InTr^e-'^^'^^ = - 5^1n(l - e'^^^O = -Trln {l - e-'^^^"^)} , (C.2) 

j 

where the Xj are the eigenvalues of the matrix K — M. In order to make sense of Zj^j we 
require the eigenvalues of i^' — M to be all positive. 

We want to find an upper bound to In Zm — IuZq. Our main result here is the following. 
(Remark: Notice that there is no (3 in ()C.3|) ). 

LEMMA C.l. Assume that < M < eK with < £ < 1. Then 

In^Af-lnZo < / Ti {{K-sM)-^M)ds = Tr ln(K)-Tr In(ir-M) < 
Jo 

(C.3) 

Proof. We write the difference on the left side of fl(I3jl as 

(C.4) 

since (e^ — 1)^^ < x^^ for x > and since M is positive semidefinite. The estimate on the 
right side of ()C.3|) follows by substituting M < sK in the denominator of the last expression 
in ()C.4|) and using the fact that is matrix monotone for x > 0. Finally, doing the s 
integral we obtain the inequality in ()C.3|) . □ 



ln(l - e) 
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We now apply this lemma to the operator in ()H.19j) . The matrix M > is given by 



Ca /27r 



(2vr) 



-r] mxmk\'/'\kf^'G{k-k',X) 



6''^ + 



|A;|2|A;'|2 



(C.5) 

with C4 = 64:aCiC/3V3 and where G{k,X) = j^,e^y^G{y,X)dy with G{y,X) given in 
()H.H2|) . (Here, X_ merely plays the role of a parameter). First, we note that, as matrices, 
M < N where 



Ca ..27r 



= ^5''{^?mx{k')H''\kr''G{k - k',x) 



TT 



(C.6) 



The requirement that < eK, as a matrix, is equivalent to the requirement that K ^/"^NK < 
el. Hence, we need to show that 



C4 .27r^3. 

TT^^ L 



(C.7) 



(C. 



j2Gik-k\x)mfm{k')f{k')<sj2\fik)\' , 

k,k' k 

for all functions f{k). The inequality 

r< o~- _ 

^(X)' E - k',X)f{k)f{k') < e J2 \fik)\' , 

k,k' k 

would clearly imply ()(I7j) . and this is implied by 

G{y,X)\J2e-''-'f{k)\'dy 

k 

<^supG(y,X)(^f / \J2^"'''fik)\'dy = 8nG,snpG{y,X)J2\fik)\\ (C.9) 



C4 27r 3 



k fc 

Thus, the condition 87rC4 sup„ x G{y, X) < e guarantees that M < eK. If we take 



e = 647rC4 = mGnaCiC/sV^ 



(C.IO) 



and use the fact that G{y,X_) < 8, the condition is satisfied. 

It remains to apply the lemma above to this particular choice of M, which (recalling 
fl3.13|) ) yields the bound 

ln(l - e) 



In Zm — In < SttG 




{^-frJ2^A{k? [ G{y,X)dy 
XAikfdk [ G{y,X)dy 
A' / Giy,X)dy 



— \\n{l-e)\K^~'^'N 



(C.ll) 

(C.12) 
(C.13) 
(C.14) 



as used in (|3.3ip . 
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